Analysis of Rotation Curves in the framework of the Gravitational Suppression model 
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We present an analysis of suitable rotation curves (RCs) of eight galaxies, aimed at checking the 
consistency and universality of the gravitational suppression (GraS) hypothesis, a phenomenological 
model for a new interaction between dark matter and baryons. Motivated by the puzzle of the core 
versus cusp distribution of dark matter in the center of halos, this hypothesis claims to reconcile the 
predictions from N-body A cold dark matter simulations with kinematic observations. The GraS 
model improves the kinematic fitting residuals, but the mass parameters are unphysical and put the 
theory in difficulty. 



PACS numbers: 95.35.+d, 04.50.+h, 98.62.Gq 
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The gravitational suppression hypothesis [ij is a phe- 
nomenological model that addresses the complex under- 
standing of the dark matter distribution on small, sub- 
galactic scales. High-resolution radio observations from 
spiral galaxies, along with their optical rotation curves 
(RCs), suggest that the dark matter is distributed in 
spherical halos with nearly constant density cores (see, 
e-g-, [1, S 01 and references therein). On the other 
hand, theoretical predictions from the well-known N- 
body ACDM simulations (e.g., [^) present a steep den- 
sity distribution profile in the centre of the halos: 
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rg is a scale radius and ps its characteristic density, in 
principle independent, but found related within a reason- 
able scatter through the halo mass, by the Bullock et al. 
equation: c = Rvir/rs ~ ^^{ i^^Mq )~°'"^^' where c is a 
concentration parameter and Rvir and M^r are the virial 
radius and mass. Mass models with a Navarro- Frank- 
White (NFW) density profile, given in Eq.(l), have two 
serious kinds of difficulty in reproducing the observed 
RCs: a) the fit is not satisfactory, i.e., Xred ^ ^ (^^^' 
e.g., 0] and references therein); b) the values of the pa- 
rameters of the best-fit mass models are clearly unphys- 
ical. In detail, the values for the halo mass result much 
higher than those we obtain from weak lensing halo mod- 
els 7jand from the analysis of galaxy baryonic mass func- 
tion [i]: Mhaio « 3 X IQI^Mq (Lb/1O"L0)1/2. In the 
same way the values of the disk mass-to-light ratio result 
much lower than those derived from colours of spirals 
BiHEl: \og {Md/Mq) « -1.6 + 1.21og(LB/L0), 
i.e., 0.7 < Md/Lb < 4. 

Several solutions have been proposed for the above is- 
sue, most of them related either to a better comprehen- 
sion of structure formation (e.g., [l^) or to new funda- 
mental physics (e.g., [l3|). Alternatively, the presence 
of noncircular motions in galaxies has been advocated to 
reconcile (up tp 70% in the Low Surface Brightness of) 
the observed kinematics with the cuspy density profile 
(e.g., [illli, but see also [B]). 

The original proposal by Piazza and Marinoni (PM) 



GraS model, instead, modifies the usual Newtonian po- 
tential of the dark matter felt by baryonic test particles 
in such a way that the NFW kinematics and the observed 
one become in agreement. According to PM, the NFW 
profile is used because GraS does not affect the dark mat- 
ter dynamics, but only the dynamics in the mixed sec- 
tor dark matter-baryons, so both primordial dark mat- 
ter perturbations and halo formation are unaffected, and 
well-known N-body simulation results can be assumed. 
The idea is adding a Yukawa contribution to the gravi- 
tational potential 
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from a hypothetical short-range interaction just between 
dark and luminous matter 
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where A is a scale range parameter. The effect is damping 
the gravitational interaction on small scales. The final 
potential is then 



4'halo — 4>Newton + OL 4>Yukawa- 



(4) 



a is a strength parameter and taken to be —1 in order 
to have the maximum possible gravitational suppression 
17| . The circular velocity is related to the potential by 

^halo — ^halo, Newton + ^halo, Yukawa = \d(l)halo / dr\. 

In PM model, for a (small) sample of RCs of Low Sur- 
face Brightness galaxies GraS was able to eliminate the 
above core versus cusp discrepancy. However, in order 
to allow a simple analytic calculation, they have taken 
a number of assumptions and approximations. In de- 
tail, the contribution to the gravitational potential from 
baryons (stars and HI disk) was neglected and the dark 
matter distribution was modeled with the simple form 
Phaio{x) = pox~^, rather than by Eq. (1). Further sup- 



port to GraS was given in [18| where the dispersion ve- 
locity of two spheroidal dwarfs (Fornax and Draco) were 
studied in this scenario. However, both large errors in 
the kinematic measurements and large geometric and or- 
bital uncertainties of the employed mass model, limited 
the relevance of their findings. 
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TABLE I: Parameters of the mass models. 



Galaxy 


Lb (Lq) 


Mass model 


Md/Lb 


71 '/■ f IK f \ 
A-lhalo (Mq) 


•2 

Xred 


rs (kpc) 


Ps (10 pcrit.o) 


Md \Mq ) 


c 


Positive results 


















ESQ 116-G12 


4.6 X lO'' 


NEW 

NFW+GraS 


0.1 

0.3 


3.8x10 
1.5x 10^^ 


2.8 

1 


14.5 ± 14 
5.1 ± 2.3 


4.0 + 6.6 
26 + 25 


(4.2 zb 11) X lU 

(1 + 1.7) X 10 


13 
26 


UGC 10981 


1.2 X 10" 


NEW 
NEW+GraS 


1.5 
0.4 


2.6x 10 
7.7x10" 


4.2 

2.5 


8 ± 2.9 
4.2 + 0.3 


13 + 9 
180 + 40 


(1.8 ± 0.3) X 10 
(4.9 + 4.4) X 10^° 


21 
55 


Negative results: unphysical parameters 














DDO 47 


10* 


NEW 


< 0.2 


7.4 X 10^^ 


1.9 


176 ± 10 


0.12 + 0.1 


< 2.3 X lO'^ 


2.8 


NEW+GraS 


0.5 


D 1 loll 

o.l X 10 


0.4 


26 + 18 


1.8 ± 1.4 


(4.5 + 2.2) X 10 


9.2 


NGC 6822 


1.6 X 10** 


NEW 


< 0.04 


1.7x10^^ 


2.3 


87 + 49 


0.19 + 0.12 


< 6.7 X 10^ 


3.5 


NEW+GraS 


< 0.02 


2.5X10 


0.5 


2.9+0.1 


24+0.7 


< 2.9 X 10 


26 


ESQ 79-G14 


2 X 10^° 


NEW 
NEW+GraS 


0.3 

0.3 


3.9 X 10 

1.1 X 10^^ 


5 

2 


o on 1 1/1 r\pi 

330 + 1400 
22.9 + 6 


0.1 + 0.49 
3.2 + 1.4 


(D.4 + 1.9) X 10 
(6 + 0.9) X 10^ 


2.6 
11.2 


UGC 8017 


4 X 10^° 


NEW 
NEW+GraS 


1 

1.1 


4.4 X lO^'' 

1.5 X 10^"* 


4 

1.6 


379 + 3600 
22 + 9 


150 + 60 
250 + 50 


(3.8 + 0.8) X 10^° 
(4.4 + 0.3) X 10^° 


51 
62 


UGC 11455 


4.5 X 10^° 


NEW 
NEW+GraS 


1.4 
< 0.2 


3.6x10^^ 
3.2x10^^ 


7.2 

3.9 


121 + 13 
13.7 + 0.5 


0.9 + 0.1 
28 + 2.6 


(7 + 2) X 10^° 
< 10"' 


7 

27 


Negative result: no change 
















M 31 


2 X 10^° 


NEW 
NEW+GraS 


6.5 
7 


1.4x10^2 
1.4x10^2 


2 

2.2 


28.5 + 1 
31 + 1.1 


2.2 + 0.1 
1.8 + 0.1 


(1.3 + 0.1) X 10" 
(1.4 + 0.1) X 10" 


10 
9.2 



In the present analysis of GraS we abandon the above 
approximations and test a wider and fairer sample of spi- 
rals. An in-depth review of the GraS model is beyond the 
scope of this Letter. Our goal is to perform a check of 
GraS. First, we assume the exact NFW profile. Second, 
we consider the baryonic contribution, so that the total 
potential is 



^■model 



'^halo 



Odisk 



^gas : 



(5) 



where the sum of the last two terms is 4>baryons- This 
leads to Vl^^^i = V^^i^ + V^^^ + V^^^. Finally, we use a 
sample of high-resolution RCs of Low and High Surface 
Brightness galaxies, in order to investigate the consis- 
tency and universality of the model. 

Our sample represents the best available RCs to study 
the mass distribution of dark matter and it has been 
used in works concerning the core versus cusp discrep- 
ancy controversy [3, [l^. The sample includes nearby 
Low and High Surface Brightness galaxies, all poorly fit- 
ted by mass models with NFW halos that also have un- 
physical values for their best-fit massparameters: DDO 
47 0; ESO 116-G12, ESO 79-G14 3: NGC 6822 [l^; 
UGC 8017, UGC 10981, UGC 11455 [2lj; M 31 0. Let 
us notice that in some cases Hq, and HI RCs are both 
available and they agree well where they coexist. More- 
over the RCs we analyse are smooth, symmetric and ex- 
tended to large radii. 

We decompose the total circular velocity into stellar, 
gaseous and halo contributions, according to Eqs.(l)-(5), 
where the latter contains the additional dark matter- 
baryons interaction. Available photometry shows that 
the stars in our sample of galaxies are distributed in 
a thin disk, with exponential surface density profile 
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FIG. 1: Galaxies in which GraS eliminates the core versus 
cusp discrepancy controversy. Y axis is the velocity in km/s. 
The solid line represents the best-fit mass model, the long- 
dashed line is the contribution of the dark matter halo, and 
the dotted and short-dashed lines are those of the stellar and 
gaseous disks. Below the RCs, we plot the residuals (Vobs — 



T.D{r) = (MD/27ri?|))e-''/^°, where is the disk mass 
and is the scale length. The circular velocity contri- 
bution is given by V|,fc(r) = {GMd/2Rd) x^B{x/2), 
where x = t/Rd and G is the gravitational constant. 
The quantity B = IqKq — Iili is a combination of Bessel 
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FIG. 2: Galaxies in which GraS does not solve the core versus cusp discrepancy controversy. The fitting values of the mass- 
to-light ratio (NGC 6822, ESO 79-G14, UGC 11455) and halo mass (DDO 47, UGC 8017) result unphysical. See Fig. 1 and 
Table I for details. 



functions [24|. The contribution of the gaseous disk is 
directly derived from the HI surface density distribution. 

In a first step, the RCs are best-fitted with the fol- 
lowing free parameters: disk mass, NFW scale radius and 
characteristic density, and scale range of GraS. Then we 
redo the analysis fixing the GraS scale range parameter 
at the mean value found of A = 3.1 kpc. Notice that the 
published mean value of PM for A is quite different from 
ours as an effect of their simplifications: A = 1.1 kpc. Our 
value is the most favourable for the PM model: different 
values of A leads to worse performance. 

The test goes against the GraS model. For the RCs 
of our sample the NFW mass halo model fails to repro- 
duce data according to the usual pattern explained in the 



introduction. Data, not surprisingly, points to dark mat- 
ter halos having inner density cores. Applying a Yukawa 
potential to the cuspy NFW halo does not solve this dis- 
crepancy. The cusp is erased and RCs are fitted very well, 
but this success is illusory in that the corresponding val- 
ues of the parameters of the best-fit mass model remain 
unphysical. In table I we show the results of the test. We 
give: the values of the parameters of the mass model and 
global properties of the galaxies. Xred calculated with 
average typical velocity errors. In bold, unphysical val- 
ues for halo mass and mass-to-light ratio, and Xred > 2.5. 
The critical density of the Universe today is taken to be 
Pcrit,o = 10~^^g/cm^. 
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FIG. 3: NFW+GraS mass model. Left: with typical values 
for the halo mass and the mass-to- light ratio. Right: with 
the Bullock et al. relation. See Figs. 1 and 2 and Table I for 
details and comparison. 

In detail, in the cases of ESO 116-G12 and UGC 10891, 
we have that GraS fits sufficiently well the RCs unlike the 
NFW, confirming that this model could work in some 
objects (see Fig. 1, Table I). 

However, in the other cases, although the fits are sat- 
isfactory, the best-fit values of the halo mass and mass- 
to-light ratio are unphysical. In fact, we expect (see 
above) the mass-to-light ratios for NGC 6822, ESO 79- 
G14, UGC 11455, to be equal to (1, 2.6, 3.5), while we 
found much smaller best-fit values (<0.02, 0.3, <0.2). In 
the same way, we expect halo masses for DDO 47 and 
UGC 8017 to be equal to (9 x 10^° Mq, 1.9 x lO^^ ^j^)^ 
while we found much bigger best-fit values (8.1 x IO^^A/q, 
1.5 X 10^'^ Mq). Furthermore, in M 31 the GraS modifi- 
cation is negligible and irrelevant (see Fig. 2, Table I). 

Let us notice that by constraining the values for the 
mass parameters within physically acceptable values, we 
obtain unacceptable fits for the GraS mass model, sim- 
ilar to those of the Newtonian NFW case. As an ex- 
ample, in UGC 8017 with Muaio = 3 x 10^^^/^ ^nd 
Md/Lb = 3Mq/Lq, GraS shows an unacceptable fit 
to data (see Fig. 3). More in general, we realize that 
for all six objects, all values of ps and Vs within their la 
uncertainties imply unphysical halo masses and/or mass- 
to- light ratios. 

We now implement the Bullock et al. concentration vs 
halo mass relation, that eliminates one parameter in the 
original NFW profile. With this relation built in, GraS 
performs even worse than before. See in Fig. 3 the case 
for DDO 47. 

In conclusion, the GraS-PM model fails to rescue the 
NFW profiles in a number of high quality well-suited 
RCs. Moreover, let us point out that there is not a pat- 
tern of this inability, so that it is presently difficult to 



understand how to modify it in order to reach its original 
goal. Then the GraS model is a rather strong hypothesis 
that does not seem solve the core versus cusp discrep- 
ancy problem of the mass distribution of the center of 
dark matter halos. 

Finally, let us remark that also in this Letter it has 
emerged that the available kinematics of galaxies is very 
constraining for non-Newtonian theories of gravity. 
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discussions. We appreciate the detailed comments from 
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version of this Letter. This research was supported by 
CAPES-Brasil (C.F.M). 
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